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1 • In trcdno tlon * In the ol&Bsical linear theory of elasticity there 1® v '* 
a widely used integral theorem called the reciprocal theorem of Betti ^ 

and Rayleigh* Sokolnikoff (ll and Love [S3 have amply illustrated the’ -• 

versatility of this theorem for elaetostatio problems while Payton J3S* 

[^) e end Boltin SsJ havs employed a dynamic version to study elastodynaaid 
problems involving moving point and line loadings c Tung [6| has 
generalised the theorem to cover the case of a linear vieooelastio ■ * ; 
eelld and his book contains references to versions of this theorem 
useful to thermoelasticity And shell theories* 

In this paper we Intend to establish a dynamic reciprocal theorem 
for a linearised theory of interaotlng media postulated in a paper by 1 - ■ 
Steel [?] • The constituents of the mixture are a linear elastio solid / 
and a linearly visoous fluid* In addition to Steel's field equations 
we use boundary conditions and inequalities on the material constants 
that' have been shewn by Atkin 0 Chadwick and Steel (8j to be sufficient 
to guarantee uniqueness of solution to initial-boundary value problems* 

The elements of the theory ©re given in seotlon 2 and two different 
boundary value problems aro considered* The reoiprooal theorem ie derived 

A 

in section 3 with the aid of the Laplace transform and the divergence 
theorem and this section Is concluded with a discussion of the special 
oases vhloh arise when one of the constituents of the mixture is absents 
As an illustration of the theorem we obtain the response of the • 
mixture occupying an infinite region and subjected to an impulsively 
applied moving point load acting on the solid constituent* The 
displacement of the solid component and the velocity of the fluid 
constituent are found and discussed* This Is the content of section 4* 
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2* Field equations for the mixfrura o W® formulate the field equations 
appropriate to & mixture of linear elastic solid and linearly viscous : 
fluid using the field equations end boundary conditions given in [?J V 
and (8]« All equations are given referred to a cartesian coordinate ; 

system x:« (x^XgpX^)? and time t® The mixture is assumed to ocoupy i 

a regular region of three-dimensional Euclidean cpac © 9 D 0 with 
bounding surface 9 S 0 The conventional Indlol&l subscript notation * 
is used to specify wctor or tensor components with an index 
appearing twice indicating a sum over 1*2 0 3« Subscripts preceded by 
a comma indicate spatial differentiation with respect to that variable, 
while time derivatives are Indicated by a dot® 

According to [7| and j8j p t ha field equations omelet *M 
following t 
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To oomploie the formulation we add to £he above th« initial andboundary 
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conditions© Thus when i*0 we require ‘ 
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w i (s,0) c w[ 2 ^b)s J?(x,0) “ .©• i*i,2®3 

for all points s in 23* while on. the boundary S we prescribe for i£0 

« t ± . 

V v i * r i 5 1 * 1,2,3 

whore are the components of the unit outward normal to S. 

Quantities appearing in t2.1) to (2.5) which are assooiated with 
the solid component of the mixture are g^ , g^ , w^ , e^ , and f^ • 
Here lo the density at time t and place x 9 pj>0 its initial value f 
w^ the displacement componcnte 9 f^ the body force components * and 
*ij 9 ^ij 1 re3 P e,otlv ® 1 y» tho strain and partial stress tensor 
, component. In the fluid p ij is the current density minus its Initial 
value | p^O, v^ the fluid velocity components, the 
rat® -of -deformation tensor, the fluid partial stress tensor, and 
the fluid body force components. The vector components in (2«5) 
are those of the diffusive resistance vector. The material constants 
# ®2 * 0% ■* P& • 9 & 9 A * y^ o y 2 and ^ are assumed to obey the 

Inequalities given in t 8] as well as the ©qualities 
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2a place of (3*6) and (2»?) we require that whan t * 0, , 

^(*,0) n V^O) « 3‘ V, 

®Wp^fe)t» B(2 SjC) ® 0* i** 9 i p 2 ? 3« ' ; 

&t all points rz of ond wh©n 25 ie en 8 # 
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Equations (2«l) 9 to (2<>7)® differ froa (3»l) to (2*7) only In 
allowing different body forces* initial conditions and surface 
conditions e notational changes are obvious and arc used for 

th© eaka of clarity in what follows o 
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by (2.1) to (2*7) and th® solutions of problem 2. given by (2«i) 9 to ■ 

' (2.7) 9 * This eolation is established with th® aid of th® LfepXao® 

■. tr&nafom and th® divergence th^oreaia 

«■ 3» Jfae Ml r aoln r oqft\ ll .t|be l 9 ^ jWtte M® begin by defining th« Laplao® fcraunefoaw 
vith respect to tins ©f & fbnotion £ (t) to bo *' * ' 
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.•• and by recalling that th® invsrs© of th® product of ^ ' 

. ■ * * 

givan by 

^{^(ajf^s) } t^(Ut)t z (x)dt . <3»2) 
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The boundary conditions (2.7) whon transferred by (3*1) bedeo* 
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New consider th® solution to (2.1) • to (2.?)( to bo given by 
; Vl^(x*t) and V^(x 8 t). Apply (3.1) to this collation* ssultiply (3«4)^ 

by rt^(x»s) and (3.4) g by ^(x 8 e), then am on 1* Integrate both 
equations ever 0 and add. Wo then ha vs 
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And this form In turn* upon using the dive-gone® 
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Equation (3*13) is the statement of the reciprocal in the transformed 
variables « & diroot inversion to real tin® yields th® final fora 
of the theorem and this can be accomplished by means of th® convolution 
(3«2)« We give several versions that ar« useful® 

i u m . 

A© .Zero initial data ® If the initial conditions (2.6} 9 (2 e 6 ) 9 car® 
homogeneous then by (3*2) and <3©13) to (3®i5) we obtain 
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B. Infinite region « If In plaoe of (2«7) and (2.7)* 

condition that velocities end stresses vanish as dietarioe increases 

t*m the origin then froa ( 3 * 16 ) there remains y : . 
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C. Single constituent * If one of the constituents is absent then * ,• •; 

froa ( 3 * 13 ) to ( 3 * 15 ) we obtain a reciprocity relation valid for 
a linear elastic solid or a linearly viscous fluid alone* l 

Suppose first that the solid is absent* Then * 0* *0 ,• • 

and the fluid equations are obtained froa (2*1) to (2*7) by : 

equating to sero the constants * 
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fluid pressure in the rest states 

Making these adjustment and redoing the boundary conditions 
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■ ttotien of a wiagturo of & linear elastic solid and visooua fluid 

As a prel&adnaxy problem w® seek w^(x*t) 0 i •« l e 2 *3 satisfying 
(2®i) to (2.5) find homogeneous Initial oonditions* 

WjfegO) s 5 & (3 e 9 0) « v^(x r 0)»® G f 1 ® i e 2 P 3 
?$(x p 0) « 0 , . (4c 

fop 'the infinite region dofi n&A by «e»<x^ $ x 2 $ x^c +00 and tSO* v 
In pl&c© of (2*7) w® require Wj(x»t) 9 v^Cxet)^ 0^ ftnAir^ k® , ; 
vanish aa (s^)^ increases without bound. • ’ ' ; 

In p&rtioul&r re consider body forces to be given by 
f(x»» 0 »t) e aj * 2; 10 ) $(* 2 - x 2Q ) - x^ 0 ) 0(t) 

g(x 6 3t @> t) « ^ £0^ - x 10 ) 6 (x 2 - x 20 ) 0(*3 “ x^q) 0(t) 

where io the unit ve©tor in th© direction and d ie the 
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Tho physical problem do scribed above corresponds to that of 
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• finding tffatXQpt) at plao© sc and at time t dua to a unit foroe 

, applied at Xq in tho direction parallel to th© axis at t ' * 0® 

The vo 6 tor a sought olearfy play tho role of Greens funotiona for 
’the theory used her©* 
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fore© loading {&©£) with g « 0 and with several restrictions on 

• th© material cohs fonts. In |9 J« the solution was given for W*i» * • 
th© diff^aslvQ resistance parameter* soro for the oases when- th© • 
fluid is iaff&saic I Q& viscous* In (10] 9 th® same problem was . '•' •• 
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Th« wave speeds o^s-Q are asseoi&ted with the elastic component 
and are defined by of - kj ^ v* * pj ^ . Kj) 1§ the 8ph8plc<a 
distance Roasured frezs the point x Q ® ‘ 

Expressions for the fluid volooity components were ale© found 

in JiO] tot sir.eo m do act intend to use then here they will net 
be repreduoedo 


I/®i ua now consider the sa&s problem with f* = 0 in ( 4 0 2 )© 
Following tho methods presented in flOj we first translate th® 
origin t® th® point s Q . then using the Fourier exponential 
transfora on eaoh of tho s'paoe variables and the Uplao® transform 
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Fourier inversion of (4»10) 1® accomplished in two stops* * 

Th© denominator of (4«10) is factored into quadra ties inAj and 
(4*10) is expanded fcy partial fractions involving these factors, 
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Inversion with rospsot to A ^ is then easily performed 0 Tfc© Inversion ■*. 
. s*«ap®ot to Ag fi A^ is next and is made easier if one exploits 

. rotational symmetry of tho expressions* 

Tho not result of those two operations leaves the function 
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. W© shall now cone! der the solid 
field produced by a moving-point 
suddenly applied ©n the solid 
t 15 0 and maintained at a eo? 
that is s m lot in (3el4) 
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Th® ' displacement field of the solid component is to bo found 
first® Since wo want to have th© displacement field and not ths 
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velocity field of th® solid component, w© go back to th® reciprocal 
statement in tha transformed variables (3.13) • ( 3 , 15 ) instead of 

utilising (3*i?)» With the aid of (4<>2) with g ® 0 and (4.17). 

w® got a direct limralen of (3®13) « (3*i5) to real tlm® in tho 
final form 
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Bventhough m are dealing with tho mixturo, th® relation* (4il8)* 
between th© displacement fields subjected to (4.2) with g rO ; 
and (4.1?) appears to bo that of th© single constituent Q.20). *♦ 
to determine th© solid displacement subject to (4.17), we 
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being understood in (4.6) and <4.7 ). 

The velocity field of the fluid component may be easily 
found by the rooipreoal relation (3.17). Considering the initial 
and regular conditions, (4.2) with 7 ■ 0. and (4.1?), we get from 
(3.17) 
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where — is the derivative with respeot to time variable 1 
frosa tho displacement (4«i5) c 

To determine V & subject to (4.17), wo substitute (4.2) with 
**® ' 

^ « 0 and (4®2J) into (4 o 20) p then wa get hgr psrfoming the . 
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BUdlttly, by employing g<x 9 * 0 *t) « ‘ / 
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^fhere a H^x-ae^fc) &ro defined to be tom ( 4 . 15 ) .*4 V’. 
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Po? the integration of (4.19) and (4.21), a careful consideration 

tauet bo given to tho behavior of the function • 

■ Sj(g) => t = | „ Sill „ g 2 (|) . « - i _ , 

1 "8 *• 1 
because the integral of (4.19) and (4.21) depend upon the aeroea 

62 ( 1 ) ®»d upon tho interval where g^(g), ff 2 (I) take 
positive values. Since the behavior of a similar function for tfes 
elastic solid was illustrated in [3S* w® omit the duplication. 

Sho, final solid displacement and fluid velocity fields are 
;,fcund to fee in polar cylindrical coordinates (r,0 9 #)#..; 
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Above wa need tho symbols R*(r 2, **j) 2 , R 0 «f(xj^s»t) 2 +(l» * .' 
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where the symbol® used are ..defined bgr 
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Stopla observation shews that the, solid displacement end fluid veld 
veleolty fiolds for the mixture exhibit components that depend upon 
the solid wave velooitles V v s and a diffusive component depending 

upon the fluid viscosities ^ «®d * 1 

Moreover, if the velocity of the moving force ie greater than ; H 

the wave velocity v, bat less than then there is a region Whoee , : 

- 1 

points B&tis^y 


R» v B t and t - - g(-^g - 1)2>0 * 


Inside this region wo hfcvo S(j/ 8 ) “ 1 and out sido this region 
S(v e ) ® 0« Therefore, th© solid displacement and fluid velocity 
t fields hav® & oomnon propagating conical wave fronts 

t m ; LtJL. „ 1 ) es o« besides the spherical one, R=v 0 t« Similarly 
* v v v£ 

if the velocity of the moving force is greater than the wave 
velooity o^, then there are two oonioal regions in which S(o 1 ) * $ 
or S(v fl ) ■ i but outside the regions 8 ( 0 ^) “ 0 and S(v B ) « 0, 
therefore 9 the solid displacement and fluid velocity fields have 
two common, propagating oonioal wave fronts besides the spherical 

one* 

The solid displacement and the fluid velocity fields have 

' ’singularities and become unbounded when R Q a 0 # «<> th * singularities 

occur atx 3 **/t 9 r *0 if V<V Q and at th® oonioal surfaces . 

‘ . x i . ; 

’ t £[(«) 2 o l] 2 a 0 if where 0 ploy the role of 

' ' 9 a Oj OP 0 8 V s • 

mally wo observe that the fluid velocity field ore of 
. Ofder0 f W both wave and diffusive somponerits. If er were aero, 
the fluid rejoice would be Identically sere. On the other hand, 
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